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ABSTRACT 

Metallic surface roughness in a nominally smooth surface is a potential indication of material degradation or damage. 
When the surface is coated or covered with an opaque dielectric material, such as paint or insulation, then inspecting for 
surface changes becomes almost impossible. Terahertz NDE is a method capable of penetrating the coating and 
inspecting the metallic surface. The terahertz frequency regime is between 100 GHz and 10 THz and has a free space 
wavelength of 300 micrometers at 1 THz. Pulsed terahertz radiation, can be generated and detected using optical 
excitation of biased semiconductors with femtosecond laser pulses. The resulting time domain signal is 320 picoseconds 
in duration. In this application, samples are inspected with a commercial terahertz NDE system that scans the sample 
and generates a set of time-domain signals that are a function of the backscatter from the metallic surface. Post 
processing is then performed in the time and frequency domains to generate C-scan type images that show scattering 
effects due to surface non-uniformity. 
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1. INTRODUCTION 

Roughness is a measure of the relative elevations of surface irregularities that result from various sources such as 
machining processes, surface damage, and corrosion. In machining applications, a machinist will compare his work 
piece to a representative sample or standard by running his finger over the work piece and sample. The drawback of this 
method is that it is qualitative with different machinists arriving at different conclusions for the same work piece. A 
contact surface-tracing machine is a quantitative method used to measure roughness that involves pulling a stylus across 
the surface. Non-contact optical methods can measure the scattering of radiation by the surface being characterized and 
yield statistical surface parameters [1,2]. One instrument for non-contact surface roughness measurement is a Total 
Integrated Scattering (TIS) instrument. This instrument measures the specular and diffusely scattered laser light from a 
surface and relates these measurements to the root-mean- squared (rms) surface roughness. This general relationship is 
[ 2 ] 


TIS = ^-= e~ i4 ^ af (i) 

where R s is the specular reflectance, R 0 is the total of the specular reflection (R s ) and the diffuse reflectance (R d ), X is 
the wavelength of light, and 8 is the rms surface roughness. In this relation, it is assumed that the light is incident 
normal to the surface. This non-contact method requires direct view of the surface being interrogated. If the surface is 
covered by paint, foam or other dielectric nontransparent material then direct interrogation is not possible. Terahertz 
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(THz) nondestructive evaluation (NDE) can potentially evaluate subsurface metallic materials because THz radiation 
can penetrate a dielectric material and inspect an underlying metallic surface. THz NDE is relatively new inspection 
technology [3, 4] that offers a non-contact and high-resolution inspection method. Its frequency regime is between 100 
GHz and 10 THz and it has a free space wavelength of 300 jam at a frequency of 1.0 THz. This radiation, like 
microwaves, is able to penetrate and inspect nonconductive material for defects and is reflected from metallic materials, 
but pulsed THz NDE returns an electromagnetic time response signal that can be interpreted in a manner comparable to 
conventional ultrasonic time-domain signals as opposed to the continuous wave methods used in microwave NDE. 
Some other uses of THz have focused on spectroscopic applications such as inspecting the chemical content of packages 
[5], inspection of artwork and diagnostics of skin burn depth and severity [6], photonic crystal characterization [7], and 
response of THz to human skin [8]. In this paper, the response of a THz NDE system is examined as a function of 
surface roughness. 


2. SETUP AND SAMPLES 

The THz system used in this work is manufactured by Picometrix, Inc and is similar to systems designed to inspect up 
to 30 cm of sprayed-on-foam insulation on the Space Shuttle external tank [9]. For this system, THz signals are 
generated and detected, using optical excitation of biased semiconductors with 70-femto-second pulses from an 800 nm 
laser. The resulting digitized time-domain signal is composed of 2048 points and has a duration of 320 ps, which 
corresponds to a propagation path length in air of 96 mm or approximately 4 inches. The transmitter and receiver units 
were in a confocal arrangement and use a 15 cm diameter, 285 mm focal length aspherical lens to concentrate the THz 
energy on the sample. The spot size diameter at the focus is approximately 3 mm. This transceiver head is scanned in a 
raster pattern over the sample while time-domain signals are digitized and stored for later processing. In the present 
work, no modifications were made to the transceiver other then scan height adjustments. The transceiver arrangement is 
shown in figure la and a sample THz signal reflected from an aluminum plate and its spectrum is shown in figure lb. 
The beam spread of the system was measured by placing a variable aperture between the THz-transmitter and focusing 
lens while monitoring the output signal amplitude. The full angle beam spread was approximately 27° and the beam 
diameter at the focusing lens was approximately 9.5 cm, thus the system f/number is approximately 3 and the Rayleigh 
limit of resolution is ~ 1 mm at a frequency of 1 THz. 



Figure 1. (a) Transceiver schematic and (b) a sample THz signal reflected from an aluminum plate and its spectrum. 







Figure 2. The surface roughness comparator plates (a) ‘Planing or Shaping’ plate and (b) ‘Casting’ plate 


The samples used for this work were comparator plates made of electroformed nickel and are standards that a machinist 
would use to evaluate his work. Two different comparator plates were used: a ‘Planing and Shaping’ plate and a 
‘Casting’ plate. The overall dimension of each plate was 5.5 cm x 12.5 cm and each roughness area was approximately 
16 mm x 25 mm. These plates covered a range of average roughness (Ra) from 0.04 jam to 100 jam, where the average 
roughness is defined as the area between the roughness profile and its mean line, or the integral of the absolute value of 
the roughness profile height over the evaluation length [10]. The plates are shown in figure 2. The ‘Planing and 
Shaping’ plate has eight separate areas of roughness. The smoothest area has an average roughness of 0.8 jam and the 
roughest is 100 jam. The roughness areas on this plate appear as parallel rulings. For the 100 jam, 50 jam, 25 jam, and 
12.5 jam areas the ruling spacing is 3.2 mm, 1.6 mm, 0.8 mm and 0.4 mm respectively. The ‘Casting’ plate has seven 
separate areas of roughness that have a range of average roughness from 50. jam to 0.8 jam. These roughness areas 
appear as randomly rough surfaces. The average roughness parameter (Ra) has an approximate relation to the rms 
roughness (8) used in equation (1), [1]. 


Ra = 0.8S 


( 2 ) 


3. OPTICAL RAY TRACING MODEL 

To gain a better understanding of how the THz system’s optical configuration might influence roughness measurement, 
an optical ray-tracing model was built and its results were compared to the experimental measurements. In this 
simulation, the sample was modeled as a mirror with its surface modified using a function that adds surface roughness 
by symmetrically rotating a surface element about its normal. To estimate the rms roughness, it was assumed that the 
surface element was rotated about it centerline and a corresponding surface height could be calculated. A series of 
surface height values were obtained by randomly varying the angle. This series of values was then used to calculate an 
rms value for the prescribed angle and then used as the rms surface roughness estimate for the ray tracing mode 
simulations. Figure 3 shows three ray tracing results for rotations of 0°, 10°, and 30° or corresponding estimated rms 
values of 0 pm, 96 pm, and 149 pm. This figure shows that as roughness increases, ray scattering increases. However, 
figure 3a shows some rays deviate from the predicted path. This appeared to be a function of the aspherical lens surface 
defined in the ray-tracing program because ray deviation does not occur when a convex-planar lens was used in place of 
the aspherical lens. To overcome this issue a large number of rays were used to average out any ray tracing errors. To 
quantify the modeling results the number of rays at the receiver was used as a measure of received specular signal 
intensity. 
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Figure 3. Ray tracing model examples with estimated surface roughness of (a) 0 jam, (b) 96 jam, and (c) 149 jam 



Figure 4. The (a) amplitude and (b) frequency images for Plate 1: ‘Planing and Shaping’ plate and Plate 2: ‘Casting’ 
plate. 


4. MEASUREMENTS AND RESULTS 

The samples were inspected with the THz system and the scan data was stored for post processing. To visualize the 
data, amplitude and frequency images were generated. The amplitude image was generated by plotting the maximum 
amplitude of each signal as a function of location and the frequency image was generated by calculating the power 
spectrum of each signal in the THz scan data and then taking the mean amplitude of the spectrum about 0.4 THz. The 
image at this frequency exhibited the best contrast between the different roughness areas. Figure 4 shows the roughness 
areas in shades of gray. The darker gray areas correspond to higher roughness and the lighter shades of gray correspond 
to less roughness. Surface roughness below 12.5 jam Ra became indistinguishable from other surfaces. 


To quantify the scan data, maximum signal amplitude and spectral content for each roughness area on the comparator 
plates were examined. The spectral content of the signals in each roughness area was examined by calculating the 
average power spectrum of all the signals in a roughness area. Figure 5 shows the average power spectrum for each 
roughness area on the ‘Planing and Shaping’ plate. The spectral differences can be identified for the 12.5 jam, 25 pan, 50 
jam, and 100 jam roughness between 0.3 and 0.6 THz. For the smoother areas, the spectra are similar and overlap. This 
illustrates that the measurement was bandwidth limited. Note that the frequency minimum at about 0.55 THz was a 
water vapor absorption line. 

TIS measurements were obtained by taking the ratio between the specular and total reflection. The specular reflection 
for each roughness area was calculated by averaging all the signals in each roughness area and then taking its maximum 
amplitude as the total specular reflection. For the total reflection, the maximum amplitude of the smoothest roughness 
area was used, in practice a smooth standard maybe used. The standard deviation of this average per roughness area was 
used to estimate a ± 0.03 TIS error. A plot of this TIS value as a function of roughness for each of the comparator plates 
is shown in figure 6 along with calculated TIS values from equation (1) and TIS values calculated from the ray tracing 
model simulations. In the model simulations and calculations, a frequency of 0.4 THz or wavelength of 0.75 mm in free 
space was used since this was approximately the center of the band-limited measurement. 

The measurements, calculations, and modeling results in figure 6 show that as roughness increases the TIS value 
decreases and that the results diverge as roughness increases. The divergence may in part be due to geometry difference 
of the samples, where the ‘Planing and Shaping’ plate has a linear roughness and the ‘Casting’ plate has a random 
roughness. In addition, errors may be introduced because a TIS instrument is typically used to measures the diffuse 
reflectance at a single wavelength over a hemisphere while the THz system measures the specular reflectance over a 
range of frequencies and over a limited angle. 



Figure 5. Average power spectrum for each roughness area on the ‘Planing and Shaping’ plate. 
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Figure 6. Measurement results as a function of rms roughness compared to calculated total integrated scattering values 
and ray tracing model results at a wavelength of 0.75 mm or frequency of 0.4 THz 


5. CONCLUSIONS 

Imaging of metallic surface roughness using THz NDE can be accomplished in a manner similar to conventional 
ultrasonic imaging. This paper shows that to quantify the data, methods of Total Integrated Scattering are plausible to 
adopt with the noted difference being a diffuse measuring system versus a specular measuring system, the THZ system 
and a single versus a band of wavelengths. The THz response to the rough surfaces showed that as roughness increased 
the scatter increased and thus the amplitude of the detected signal decreased. This can be seen in the images and the TIS 
versus Roughness plot where surface roughness below 12.5 qm Ra or ~10 qm rms became indistinguishable from other 
surfaces. The best correspondence between measurements, formula, and modeling appear to be in the rms range from 
~10 qm to ~ 40 qm, above this range, the results diverge. This divergence may be due to instrument, modeling, and 
sample geometry differences. With modeling and technique refinements, the useful roughness measurement range may 
be extended. 
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